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Abstract. Inversions to determine the squared isother- 
mal sound speed and density within the Sun often use 
the helium abundance Y as the second parameter. This 
requires the expHcit use of the equation of state (EOS), 
thus potentially leading to systematic errors in the results 
if the equations of state of the reference model and the 
Sun are not the same. We demonstrate how this potential 
error can be suppressed. We also show that it is possible 
to invert for the intrinsic difference in the adiabatic expo- 
nent Fi between two equations of state. When applied to 
solar data such inversion rules out the EFF equation of 
state completely, while with existing data it is difficult to 
distinguish between other equations of state. 
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1. Introduction 

Solar oscillation frequencies can be inverted to determine 
the internal structure of the Sun. This is generally done by 
relating the frequency differences between a solar model 
and the Sun to differences in the structure by linearising 
the oscillation equations under the assumption of hydro- 
static equilibrium. If, for example, we express the frequen- 
cies in terms of the squared adiabatic sound speed c^ and 
the density p, the result can be written 
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(e.g. Dziembowski et al. 1990). Here 6uji is the difference 
in the frequency Ui of the ith mode between the solar 
data and a reference model. The kernels K^2 and K^ 2 
are known functions of the reference model which relate 
the changes in frequency to the changes in c^ and p re- 
spectively. The term in F^ results from the near-surface 
errors in the model, such as the assumption of adiabatic 



oscillations; Ei is the inertia of the mode, normalised by 
the photospheric amplitude of the displacement. 

The kernels for the (c^ , p) combination can be eas- 
ily converted to kernels for others pairs of variables like 
(Fi,p), (u,ri) with no extra assumptions (cf. Cough 
1993); here u = p/ p is the squared isothermal sound speed, 
p being pressure. However, very often in addition to the 
oscillation equations it is assumed that the equation of 
state is known; this may be used to construct the kernels 
for u or p and Y . The implicit assumption that the equa- 
tions of state in the Sun and reference model are the same 
leads to potential systematic errors in the inversion for u 
or p. 



2. Formulation of the inverse problem 

The conversion of the kernels for {(? , p) to those for 
(u, Y) uses 5\nc? = SlnTi + Slnu. Had the equation of 
state been known, 6T1 could have been determined from 
Fi = Ti(j), p, {Xi}), where {A,} is the composition. In 
fact, 6T1 contains an additional term, viz. the intrinsic dif- 
ference ((5Fi/Fi)int at fixed (p, p, {Aj}) between the true 
and the model equations of state. Characterizing the com- 
position by the abundances by mass Y and Z of helium 
and heavy elements and assuming an unchanged heavy- 
element abundance Z , we therefore get 
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Note that differences in Z, or in the relative composition 
of the heavy elements, will also appear in ((5Fi/Fi)int- 

Equation (2) can now be used to rewrite Eq. (1) in 
terms of 6u/u and SY, by expressing the terms in Sp/p 
and 5p/p in terms of 5u/u by means of the equation of 
hydrostatic support. This expression contains a contribu- 
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tion from the intrinsic difference in Fi weighted by the 
kernel for c^ at constant p. Thus, the full equation is: 
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In the expression normally used to invert for {u, Y) the 
term in {5Ti/Ti)-^^^ is ignored; this clearly introduces a 
systematic error, if the equations of state are in fact dif- 
ferent. A similar argument also holds for density inver- 
sions using the pair (p, Y) . Here too the contribution from 
{5^il^i)i-a_t should be taken into account. We also note 
that using Eq. (3), one may directly invert for {5Ti/Vi)-^^^. 



3. Inversion technique 

We have used the Subtractive Optimally Localised Avera- 
ges method of Pijpers & Thompson (1992), adapted to 
inversion for structure differences (e.g. Basu et al. 1996). 
The principle of the inversion technique is to form linear 
combinations of Eqs (3) with weights Ci(ro) chosen such 
as to obtain an average of, for example, 5u/u localised 
near r = tq while suppressing the contributions from the 
remaining terms in Eqs (3), including the near-surface er- 
rors. In addition, the statistical errors in the combination 
must be constrained. 

To invert for 5u/u the coefficients Ci are chosen to min- 
imise 

jn^c^Kly^T] dr + fS^I lj2cMr)Kl.J dr 
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with the constraint that the averaging kernel be unimodu- 
lar, i.e., 
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(5) 



Here, T(ro, r) is a target averaging kernel, chosen to be a 
Gaussian of unit area centred at tq. Eij is the covariance 
matrix of errors in the data. The parameters Pi and (32 
control the contributions of 6Y and ((5ri/ri)int, respec- 
tively, and fj, is a trade-off parameter which controls the 
effect of data noise. The function w{r) is a suitably cho- 
sen, increasing function of radius, which ensures that the 
contributions from the second and third terms from the 
surface layers are suppressed properly. 

To reduce the influence of near-surface uncertainties 
we apply the additional constraints that 
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where the $a are B-Splines with a suitably scaled argu- 
ment (cf. Dappen et al. 1991). 

To carry out an inversion for {STi/Ti)ijit one minimises 
instead 

j \11<'^K^,P-A dr + PiJ r£cMr)Kly] dr 



(7) 



We have used four solar models for this work. All 
models have been constructed with OPAL opacities (Igle- 
sias, Rogers & Wilson 1992) at temperatures higher than 
10^ K and Kurucz tables (Kurucz 1991) at lower tem- 
peratures. The models have been constructed with dif- 
ferent equations of state — Livermore (OPAL) (Rogers, 
Swenson & Iglesias, 1996) MHD (e.g. Mihalas, Dappen 
& Hummer 1988), EFF (Eggleton, Faulkner & Flannery 
1993) and CEFF (cf. Christensen-Dalsgaard & Dappen 
1992). The properties of the models, identified by the EOS, 
are summarised in Table 1. Model OPAL is Model S of 
Christensen-Dalsgaard et al. (1996). 

Table 1. Solar models used; dcz is the depth of the con- 
vection zone, and Tc and pc are central temperature and 
density 



EOS 


Z/X 


dcz/R 


(106 K) 


Pc 

(g/cm^) 


OPAL 


0.0245 


0.2885 


15.67 


154.2 


MHD 


0.0245 


0.2876 


15.67 


154.5 


CEFF 


0.0248 


0.2863 


15.68 


155.0 


EFF 


0.0248 


0.2852 


15.74 


157.2 



We use solar oscillation data obtained by the LOWL 
instrument during the first year of data collection (Tom- 
czyk et al. 1995, Schou and Tomczyk, in preparation). 
The dataset consists of modes of degrees to 99 in the 
frequency range 1 to 3.5 mHz. The observed modeset and 
errors were also used in tests for solar models, in order to 
get realistic properties of the inversion. 



4. Results 

In Fig. 1 we show inversion results, without and with sup- 
pression of iSTi/Ti).^^^, for u in Model MHD with Model 
OPAL as the reference model. The resolution of the inver- 
sion is the same in both cases. Note that when the contri- 
bution from (<5ri/ri)jj^^ is not constrained, the results are 
not very accurate. The results improve dramatically, par- 
ticularly in the core, when the intrinsic difference is taken 
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Fig. 1. The inversion for the squared isothermal sound speed 
(«) difference between model MHD and model OPAL. The 
solid line is the exact difference and the points are the difference 
obtained by inverting the frequency differences between the 
models, a Inversion results when the intrinsic difference in Fi 
between the OPAL and MHD equations of state is ignored. 
b Inversion results when the intrinsic difference is taken into 
account. The vertical error-bars arc la propagated errors 




Fig. 3. Averaging kernels for the inversion for the intrinsic 
((Sri/ri)int with the LOWL mode set and the OPAL model as 
the reference model 




Fig. 2. The same as Fig. f , but for the density difference be- 
tween models MHD and OPAL 



into account. However, the price paid for increased accu- 
racy is decreased precision, as reflected in the increased 
error-bars. Indeed, it is evident that the use of the data to 
suppress the possible error in the EOS reduces the amount 
of information available for the determination of Su/u and 
hence results in larger errors if the resolution is kept ap- 
proximately the same. The results for density inversion are 
shown in Fig. 2. We note that the errors in the inversion 
corrected for a possible inconsistency in the equation of 
state are much reduced if accurate data on higher-degree 
modes are available, as is the case, e.g., for the frequen- 
cies obtained by the SOI/MDI experiment on SOHO (cf. 
Kosovichev et al. 1997). 

Although it is useful to be able to suppress the ef- 
fects of errors in the equation of state when inverting for 
u, it is evidently of greater interest to obtain a localized 
measure of these errors, i.e., to invert for the intrinsic 



{STi/Ti)int between the equations of state of the Sun and 
the model. To illustrate our ability to achieve such local- 
ization, Fig. 3 shows averaging kernels for the inversions 
for {STi/Ti)int through minimization of expression (7). 
Note that for rg > 0.5Rq the averaging kernels are quite 
well localized, indicating that reliable inversion is in fact 
possible. The averaging kernels are not as small near the 
surface as one would hope for; inclusion of higher-degree 
modes would substantially improve the behaviour in this 
region. 

Fig. 4 shows the inversion for the intrinsic differences in 
Fi , using the Models MHD, CEFF and EFF as test models 
and Model OPAL as reference. For comparison are shown 
exact differences resulting from differences in the equation 
of state, evaluated at fixed p, p, and Y in Model OPAL. 
The inversion of the frequencies clearly successfully repro- 
duces even the subtle intrinsic EOS differences between 
the MHD and OPAL formulations, although the statisti- 
cal errors are fairly substantial compared with these dif- 
ferences, at least beneath the dominant ionization zones. 

Given the success of this test on artificial data, we may 
consider differences between the solar and the model equa- 
tions of state, as obtained from analysis of the observed 
frequencies. Fig. 5 shows the resulting intrinsic differences 
in Fi between the Sun and the four models of Table 1. It 
is evident that the EFF equation of state is inconsistent 
with the data. With the current level of errors, it is dif- 
ficult to distinguish between the other three equations of 
state. Our ability to do so would be greatly improved by 
analysis of higher-degree data, since we may expect that 
the dominant differences in the equations of state are close 
to the surface of the Sun. 
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Fig. 4. The results of inversion for the intrinsic Fi difference 
between a EFF, b CEFF and c MHD models and Model 
OPAL. The solid line is the exact difference and the points 
are the difference obtained by inverting the frequency differ- 
ences between the models. Note that the scale in Panel a is 
much larger than that in Panels b and c 
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Fig. 5. The intrinsic Fi difference between the Sun and the 
EFF, CEFF, MHD and OPAL models obtained by inversion 
of LOWL Year-f data. Note the difference in scale between 
panel a and the other panels 

5. Conclusions 

We have shown that inversions for the squared isothermal 
sound speed u and the density p may suffer from systema- 
tic errors vifhen based on the common implicit assumption 
that the equations of state in the Sun and the reference 
model are the same. These errors can be removed by sup- 



pressing the contribution from the intrinsic difference in 
Fi to the frequency difference. However, this is achieved 
at the price of an increase in the propagated errors. 

We also show that we can successfully invert for the 
intrinsic difference in Fi between the currently available 
equations of state. This differs from the analysis by Elliott 
(1996) who investigated the EOS in terms of the total dif- 
ference between the solar and the model Fi. Inversions of 
solar oscillation frequencies show that the EFF equation 
of state can be ruled out by direct inversions. With the 
current level of data errors, it is difficult to judge the sig- 
nificance of the differences between the solar equation of 
state and the CEFF, MHD and OPAL equations of state. 
We hope, however, that as more precise data and data on 
high-degree modes become available, this method can be 
used as a direct test of the solar equation of state. 
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